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Taking liquid xenon detectors as a case study, the importance of a robust recoil energy calibration 
as a prerequisite to a search for light-mass Weakly Interacting Massive Particles (WIMPs) is em- 
phasized. Important shortfalls in the analysis of existing measurements of the relative scintillation 
efficiency (£ c ff) and ionization yield (Q y ) for nuclear recoils in liquid xenon are described, leading 
to the conclusion that recent attempts to extract light- WIMP sensitivity limits from the XENON 10 
and XENON100 detectors are premature and overly optimistic. 
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I: INTRODUCTION 

Direct searches for dark matter WIMPs are remark- 
ably difficult experiments where an attempt is made to 
extricate rare signals from backgrounds orders of mag- 
nitude more frequent. An additional complication arises 
from eventually having to demonstrate, if at all possible, 
that any irreducible events are more likely the result of 
WIMP interactions than a less exotic source. In the case 
of light- WIMPs below a mass of -10 GeV/c 2 , the sit- 
uation is markedly more complex, given the modest nu- 
clear recoil energies expected from their interactions (few 
keV or even sub-keV, depending on WIMP and target 
masses) , typically very close or below the energy thresh- 
old of the detectors. This is an energy region for which 
models of detector response and expected backgrounds 
can break down, background rejection becomes inade- 
quate or non-existent, and events with a non-radioactive 
origin (microphonics, electronic noise) may be taken for 
a WIMP signal. Recent times have seen a renewed inter- 
est in this WIMP mass region. The drive has been the 
detection by the DAMA/LIBRA collaboration of an an- 
nual modulation in the low-energy counting rate of their 
scintillators, in agreement with expectations from a dark 
matter galactic halo composed by light dark matter par- 
ticles [1]. A second source of interest comes from particle 
phenomenologies that generate dark matter candidates 
precisely in this mass range [2]. The situation has been 
compounded by the observation of a so far unexplained 
background excess in CoGeNT [3] , a detector targeted to 
search for light WIMPs, and more recently by an excess 
of low energy recoils in CRESST bolometers [4, 5], both 
in possible agreement with the DAMA/LIBRA results if 
the WIMP hypothesis is adopted [6]. 

The purpose of this note is to provide a critical inspec- 
tion of the methods that have been proposed or employed 
to understand the response to low-energy nuclear recoils 
and to derive a light- WIMP sensitivity from liquid xenon 
(LXe) devices. Important deficiencies are found in these 
methods. The conclusion drawn is that the knowledge 



necessary to allow a reliable exploration of this WIMP 
mass region is presently absent for LXe-based detectors. 

The outline of this paper is as follows: Section II 
deals with the difficulties in establishing a reliable en- 
ergy scale for few keV r recoils in LXe, in particular via 
the so-called best-fit Monte Carlo method. Recent en- 
ergy scales proposed using this method are compared to 
existing expectations and to independent measurements 
of ionization yield by low energy Xe + impact on surfaces, 
finding a large disagreement in both cases. Well-known 
physical processes presently neglected by this method 
are described. A commentary is provided on how their 
inclusion should affect the energy scale, easing the dis- 
agreements described. A brief attempt at the interpre- 
tation under different energy-scale scenarios of a low- 
energy feature present in a recently released XENON10 
ionization spectrum is offered. Section III describes an 
important systematic effect neglected in the analysis of 
all XENON10 and XENON100 measurements of rela- 
tive scintillation efficiency 1 (£ c ff) using monochromatic 
neutron scattering. A discussion of its impact on the 
energy dependence of £ c ff is provided. Separately, the 
low quality of the experimental data in recent attempts 
by XENON100 to characterize £ e ff is pointed out. Fi- 
nally, these concerns are illustrated with data from on- 
going quenching factor measurements at the University 
of Chicago. Section IV contains the conclusions. 

II: UNCERTAINTY IN RECOIL ENERGY SCALE 

In a recent workshop presentation [7], the XENON10 
collaboration released a preliminary spectrum of events 
in their detector as a function of ionization yield, i.e., the 



1 For the purpose of the present discussion C e ff can be considered 
the equivalent of the more broadly employed concept of "quench- 
ing factor", in this case the ratio between the yield of direct (SI) 
scintillation light produced by a nuclear recoil and that from an 
electron recoil of the same energy, at zero drift field. 
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number of electrons extracted from LXe via the applica- 
tion of an electric drift field Ed— 0.73 kV/cm. These elec- 
trons are further accelerated through the gaseous phase 
of the detector, producing a considerable VUV light emis- 
sion via electroluminescence (the so-called S2 light), re- 
sulting in the registration of ~25 photoelectrons in the 
photomultipliers per each electron extracted from the liq- 
uid phase. The high gain provided by the electrolumines- 
cence permits the detection of single electrons with good 
resolution. It is therefore expected that this ionization 
spectrum should reach down to a threshold in recoil en- 
ergy of 0(1) keV r , allowing a search for light WIMPs. 
A dark matter search based exclusively on an analysis of 
the LXe ionization spectrum could in principle circum- 
vent many recent discussions about the energy depen- 
dence of the relative scintillation efficiency £ e ff, and how 
it impacts LXe sensitivity to light WIMPs [8, 9]. While 
electron-recoil background rejection is lost when consid- 
ering just the ionization (S2 light) spectrum, the low- 
energy counting rate 2 achieved in XENON10 (<1 count 
/ keV r kg day, Fig. 1) should be sufficient to investigate 
the remaining favored DAM A/LIBRA region [6], if the 
energy threshold is indeed as low as expected. However, 
an important requirement for the analysis of such ioniza- 
tion spectra in terms of sensitivity to light WIMPs is the 
establishment of a reliable correlation between the recoil 
energy scale and the ionization yield measured. During 
their recent presentation, XENON10 collaborators em- 
ployed a new proposed recoil energy scale to generate the 
spectrum in Fig. 1 (top). If confirmed, it would amply ex- 
clude the WIMP phase-space regions presently compati- 
ble with a dark matter interpretation of DAMA/LIBRA, 
CoGeNT and recent preliminary CRESST data [6], at 
the very least for the standard present understanding of 
WIMP-nucleus interactions. 

Reservations were expressed in [7] about the trustwor- 
thiness of the recoil energy scale employed to arrive at 
this spectrum. The vast extent to which this caution is 
advised is emphasized here, if the goal is to obtain reli- 
able light- WIMP limits: Fig. 2 displays the most recent 
attempts [7, 11] by the XENON10 collaboration to es- 
tablish a correlation between the ionization yield and the 
recoil energy scales, represented by solid color lines and 
matching color bands for the claimed one-sigma uncer- 
tainties (statistical error only). The ionization yield vs. 
recoil energy is typically represented by the normalized 
quantity Q y (units of electrons/keV r ). For clarity, here 
the actual measured quantity (Q, charge yield in number 



2 This background level (~0.06 ckkd at 20 keV ee before fiducial- 
ization, ~0.01 ckkd following it), has been widely advertised in 
recent XENON100 presentations as being two orders of magni- 
tude lower than in any competing dark matter detectors. In 
reality, it is comparable to the past generation of intrinsic ger- 
manium detectors (^0.03 ckkd, [10]). 
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FIG. 1: Top: Spectrum of events in the XENON10 detector 
based on ionization yield only, as presented in [7] . The energy 
scale was defined by a "best-fit Monte Carlo" method. The 
original spectrum of events vs. electrons extracted from LXe 
can be found in [7]: the lowest bin shown here corresponds 
to an ionization yield of 3-4 electrons. The color code cor- 
responds to different fiducial volume cuts [7]. The predicted 
signal for two example light WIMPs (dotted lines) is over- 
lapped on the spectrum. Bottom: Similar to top panel, but 
using the dotted expectation curves in Fig. 2 to define the 
recoil energy scale (see text). 



of electrons) is shown on the horizontal axis. All these 
energy calibrations originate in variations around a com- 
mon method: a comparison is established at some point 
between experimental data obtained during exposure to 
a neutron source (commonly AmBe) and a Monte Carlo 
simulation of the expected spectrum of events vs. recoil 
energy. In the case of the latest attempt ([7], red line), 
this is accomplished by using a spline with mesh points 
fixed at arbitrary recoil energies (1,2,4,8. . . 128,256 keV r ) 
to describe Q y , allowing Q y to float unconstrained until 
the best possible match between the experimental ion- 
ization spectrum and the simulated spectrum of events 
vs. recoil energy is reached 3,4 . The so obtained Q y de- 
fines the correlation between recoil energy and ionization 
scales. As pointed out in [8], the original disagreement 



No information was provided in [7] on how discrepancies in the 
overall rate of simulated and measured events are dealt with 
during this procedure, or about the selection of fitting range. 
These choices can, on their own, alter the extracted energy scale. 
A similar method has been used in the past to obtain C e ff [12]. 
At least there, the overall event rate normalization between sim- 
ulation and data was treated as a free parameter, a highly ques- 
tionable approach (see Sec. Ill) able to affect the low energy C e ff 
and Q y . 
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FIG. 2: Colored lines: correlations between recoil energy scale 
and ionization yield proposed by the XENON10 collaboration 
over the last two years, using best-fit Monte Carlo methods. 
The colored bands represent the claimed one-sigma statisti- 
cal uncertainties (see text). Dotted lines correspond to the 
most recent expectations, based on the formalism presented 
in [14] (one based on "Columbia" detector data, the second 
on the "Case" detector). The dashed line indicates an earlier 
attempt to predict these [15]. All expectations are for Ed= 
0.73 kV/cm. The data points correspond to measurements by 
Manzur et al. [16, 17]. Red arrows indicate the position of the 
endpoint for the low-energy rise in the spectra of Fig. 1 (see 
text). Inset: correlation between 2.8 MeV neutron scattering 
angle and recoil energy deposited in LXe in the measurements 
by Manzur et al. [16]. 



between simulations and LXe data is typically severe at 
low recoil energy: spectra of direct (SI) scintillation dis- 
play a lack of response to AmBe neutron-induced recoils 
below few keV r when compared to the expectations, re- 
gardless of C e ff adopted [13]. Possible reasons for this are 
described below. While no mention of this initial state of 
affairs for the S2 light (ionization) spectrum was made in 
[7] , the limited information provided indicates that this is 
the case too. The XENON10 collaboration refers to these 
methods collectively under the unnerving denomination 
"best-fit Monte Carlo" [11]. 

A reader familiar with the energy calibration of radi- 
ation detectors will readily identify the risks involved in 
this approach: all uncertainties in the inputs to the simu- 
lation, limitations to its accuracy 5 , any experimental sys- 



5 To this author's knowledge, Geant4 (used in [7]) has never been 
validated or verified in the ~few keV to sub-keV recoil energy 
region, having historically exhibited difficulties in correctly mod- 
elling low-energy neutron transport. Several important improve- 
ments to the code in this respect are still in the making [18] . 



tematic effects, and most importantly, the effects of any 
physical processes not included in the simulation (such as 
the kinematic cutoff discussed below), will be automati- 
cally reflected as distortions in the inferred recoil energy 
scale. In other words, in this method, the reasons for 
any initially noticeable disagreements between simula- 
tion and experimental data are not investigated, instead 
the recoil energy scale is used as a buffer to absorb these. 
Taking as a reference a yield of 10 electrons (the approxi- 
mate endpoint for the observed low-energy rise in Fig. 1), 
the low-energy extrapolation of the colored curves in Fig. 
2 indicates that this method of calibration results in re- 
coil energy assignments for it anywhere in the ~0-2 kcV r 
interval. No explanation is provided for the origin of the 
inflexion points and non-proportionality (particularly un- 
realistic at high energy) noticeable in those curves. The 
trend of a marked increase in Q y towards decreasing re- 
coil energy in two of these curves (red, blue) is against all 
expectations: as pointed out in [8], a kinematic cut-off at 
^40 keV r is expected for nuclear recoils in LXe, below 
which ionization should be adiabatically quenched. The 
same can be said of primary SI scintillation, to the extent 
that it is expected to be mainly mediated by ionization 
[8] . Indeed, the observed traditional deficit of low-energy 
recoil events in LXe under AmBe neutron irradiations 
mentioned above is easily understood under this consid- 
eration. 

The origin for this kinematic cutoff is in simple two- 
body kinematics: when the largest possible energy im- 
parted to a valence electron by a slow-moving recoiling 
ion falls below the minimum excitation energy of the sys- 
tem (a 9.3 eV bandgap in LXe), a progressive quenching 
of scintillation and ionization should be observed (other 
secondary processes such as for instance potential elec- 
tron emission from Auger de-excitation in an ionized pro- 
jectile, whenever allowed, or transient autoionizing quasi- 
molecules can still play a role, leading to a smooth rather 
than sharp cutoff [19, 20, 26]). Contrary to an opinion 
recently expressed by the XENON100 collaboration [9], 
the basic principles behind this intrinsic limitation are 
well-known, with abundant references in the experimen- 
tal and theoretical literature on transport of slow ions 
[19, 20, 26]. This quenching has been observed at the 
predicted kinematic cutoff for hydrogen recoils in organic 
scintillator [19]. Preliminary evidence for it in Cs and I 
recoils in CsI[Na], also at the expected cutoff energy, is 
provided in Sec. Ill of this paper. Due to an unfavor- 
able combination of nuclear mass and electron band-gap 
energy, LXe-based dark matter detectors should be par- 
ticularly affected by this limitation: the cutoff for Na 
recoils in NaI[Tl] and Ge recoils in germanium crystals is 
predicted to appear below the threshold of present dark 
matter detectors [8]. Other less well-known radiation ef- 
fects left out by the best-fit Monte Carlo method could 
be listed: for instance, as recoil energies become very 
small, a diminishing recoil track length relative to the 
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range of ionized electrons may lower the chances of re- 
combination, helping any charge still being generated to 
escape and contribute to the S2 light [16]. 

It cannot be overemphasized that none of these mi- 
croscopic physical processes mediating the generation of 
information carriers (free electrons, primary scintillation 
photons) following a nuclear recoil are included in popu- 
lar simulation packages such as Geant4 or MCNP, which 
stop at generating a simple distribution of recoil ener- 
gies. Just for this reason alone, the best-fit Monte Carlo 
method should be expected to generate extraneous struc- 
ture in the correlation between recoil energy and ioniza- 
tion yield, aberrating the true relationship between these 
two magnitudes 6 . More specifically, an attempt to ex- 
pand the method by including these processes prior to 
the reconciliation of simulation and data should alter the 
behavior of the colored curves in Fig. 2: for instance, the 
adoption of an adiabatic term [19] to include the reduced 
efficiency for electronic excitation below kinematic cutoff 
should bring those curves closer to or beyond the expec- 
tations discussed next (an adiabatic term would rather 
naturally account for any observed decrease in response 
to AmBe neutron recoils at low-energy). The best- fit 
Monte Carlo method, as it stands, cannot be accepted 
as a substitute for a genuine effort to understand these 
underlaying physical processes and the way they impact 
the recoil energy scale. 

During this last attempt at defining a recoil energy 
scale, the XENON10 collaboration neglected any men- 
tion to the subject of expectations. These exist, and have 
been discussed by this collaboration and other workers 
before. The dotted lines in Fig. 2 show the expected cor- 
relation between recoil energy E rec and ionization yield, 
calculated following the relationship described in a re- 
cent XENON10 publication [14], Q{E d ) = 0.2 CE rec /W e , 
where C is the Lindhard nuclear quenching factor 7 (ex- 



6 Leaving aside the limitations of the best-fit Monte Carlo method 
in its present form, the derivation of the new WIMP exclusion 
plot presented in [7] begs justification: the one-sigma error bars 
on Q y , converted here to the red band in Fig. 2, are exclusively 
statistical. They are therefore vastly underestimated, not in- 
cluding the uncertainties in the energy scale generated by this 
method. The broad dispersion in the curves depicted in Fig. 2 
gives an idea of the magnitude of this (hard to quantify) un- 
certainty in the energy scale. Even if this one-sigma band is 
accepted at face value, the customary WIMP exclusion contour 
should be extracted from a spectrum with recoil energy scale de- 
rived from a conservative excursion from a central Q y value, as 
allowed by 90% C.L. uncertainties. Taking again the ~10 elec- 
tron endpoint as a reference and using Fig. 2 as a guide, it is 
possible to estimate that the 90% C.L. dark matter sensitivity 
claimed in [7] should be relaxed by more than an order of mag- 
nitude. To add to the confusion, a markedly more conservative 
(by a factor ~40) exclusion plot based on the XENON10 S2 spec- 
trum was recently presented in [21], using an identical Q y , i.e., 
the same proposed recoil energy scale as in [7]. 

7 A discussion of the relationship between C e ff and C (the second 



tracted here from SRIM2010 [23]), W e = 15.6 cV is the 
average energy spent by an electron recoil to form an ion- 
electron pair [22] , and the charge yield at a drift field Ed 
relative to that at infinite field Q(Ed) /Q(oo) = 0.2 was 
measured at 56.5 keV r in [14]. These expectations, based 
on the Lindhard-Scharff-Schiott theory through their de- 
pendence on C, generate a Q y roughly constant in en- 
ergy, i.e., a proportionality between recoil energy and 
ionization. It must be kept in mind that the effects of 
the kinematic cutoff may not be completely accounted 
for in such models, as noticed experimentally in [19]. In 
other words, more sophisticated models generating ex- 
pectations tending to the horizontal at few keV r in Fig. 
2 can be constructed. 

There is an evident tension between these expecta- 
tions and the ever mutable low-energy scales so far pro- 
posed by XENON10, most certainly the result of the ar- 
tificial structure that the present best-fit Monte Carlo 
method is expected to introduce. However, it should 
be kept in mind that this expectation is itself subject 
to considerable uncertainty, and therefore should be re- 
garded as simple guidance. For instance, the value of 
Q(Ed)/Q(oo) = 0.2 employed to generate the dotted 
curves in Fig. 2 may have an unknown variation over 
the broad energy range depicted in this figure. Similarly, 
while the latest version of SRIM [23] , used here to gener- 
ate C, offers a considerably better agreement with exist- 
ing LXe measurements than previous ones [24] , it should 
be kept in mind that SRIM predictions are semiempirical, 
i.e., biased by low-energy quenching factor measurements 
that may in turn be flawed. In particular, an absolute 
cutoff for £ c ff at few keV r , a possibility strongly argued 
for in Sec. Ill, would suggest an expectation curve in Fig. 
2 tending to the horizontal at this cutoff energy. 

Just to illustrate the incipiency of the knowledge about 
the processes mediating the generation of ionization and 
scintillation by low-energy recoils in LXe, a dashed line 
in Fig. 2 shows the expectations for nuclear recoils put 
forward by XENON10 as recently as in [15], which were 
inferred from observations made using alpha particles. 

Besides the colored curves in Fig. 2, generated via the 
best-fit Monte Carlo method, few other measurements of 
Q y exist. Measurements of Q y in [14], not shown in Fig. 
2 for clarity, approximately follow the blue curve over 
the range 20-100 kcV r . However, the recoil energy scale 
used to extract Q y in [14] was derived using a debat- 
able choice of £ c s. The methodology of the calibrations 
that generated this £ c ff curve will be sharply criticized 
in Sec. Ill of this paper. A rapidly decreasing £ c ff to- 
wards zero recoil energy like the alternative proposed in 
Sec. Ill would result in changes to the recoil energy scale 
that would make the Q y derived in [14] tend towards the 



sometimes referred to as q nc i) can be found in [16]. 
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dotted expectation lines rather than the blue curve in 
Fig. 2. Erroneous conclusions from one LXe calibration 
can rapidly propagate to another, such is the complex- 
ity of the relationship between quantities and physical 
processes in LXe detectors. 

The data points shown in Fig. 2 deserve special at- 
tention. They correspond to measurements by Manzur 
et al. [16] at a value of Ed similar to that used in the 
XENON10 detector, with the important difference that 
their energy scale is in principle known, independently 
determined by the kinematics of the elastic scattering of 
monochromatic 2.8 MeV neutrons from a xenon target. 
The apparent trend in those measurements is one of an 
increasing Q y with decreasing recoil energy. Not only this 
clashes with the expected behavior below kinematic cut- 
off, but it also poses an interesting number of questions: 
the low-energy extrapolation of these datapoints in Fig. 
2, if confirmed, would imply that the low-energy rise be- 
low ~10 electrons observed in the XENON10 ionization 
spectrum (Fig. 1) cannot originate in nuclear recoils. An 
alternative would be a process not the immediate result 
of a particle interaction, for instance, an "spontaneous" 
multiple electron emission. Spontaneous single-electron 
emission from LXe is known to follow large energy depo- 
sitions by gamma rays, with a half-life of ^100 fis, and 
is possibly due to photoionization [25]. Delayed emis- 
sions are not uncommon in some detecting media: for 
instance, a few inorganic scintillators (e.g., CsI[Tl]) are 
notorious for an "afterglow" from long-lived phosphores- 
cent states, involving the release of single scintillation 
photons. However, it is very hard to envision an atomic 
energy-storage mechanism in LXe that would lead to the 
delayed emission of up to ~10 electrons. The high gain 
afforded by the electroluminescence clearly points at a 
multiple rather than single electron emission. The alter- 
native left is to ascribe an origin in minimum ionizing 
particle interactions to this observed rise in rate below 
10 electrons. The value of Q(Ed)/Q(oo) measured in 
[14] for electron recoils suggests that this rise would then 
have an endpoint at an electron-equivalent energy of a 
mere ~0.25 keV ee . Compton scattering from gammas 
is not expected to generate any such features (the Klcin- 
Nishina relation does not include a highly forward-peaked 
excess). Beta decay or Bremsstrahlung interpretations 
concentrated this low in energy are not plausible. The 
chance of this signal arising from degraded surface activ- 
ity seems to be excluded by its survival under different 
fiducial volume cuts (Fig. 1, [7]). No other meaningful 
processes come readily to mind. While an understanding 
of this low-energy feature would be desirable, the appar- 
ent low-energy trend in the Q y measurements by Manzur 
et al. limits the interpretations. 

The measurements of Q y by Manzur et al., taken at 
face value, generate a third concern. The processes gov- 
erning electron emission by slow ions immediately prior 
and during impact on surfaces are not identical to those 



expected from recoils in bulk LXe, but the dominant 
mechanisms are in common. Therefore, it is educational 
to contrast the several tens of electrons generated by a 
1 keV r Xe recoil that would be implied by the extrap- 
olated low energy Q postulated by Manzur et al., with 
the observed electron yields of O(10 -3 ) e~/ion generated 
by a Xc + ion impact of the same energy on W or Au 
electrodes [26]. In this energy regime the ion range and 
the mean escape depth of electrons from the electrode 
are comparable (few nm, [26]), i.e., the measured yield 
is representative of the actual total ionization generated. 
This very small yield is of the same order of magnitude of 
what would be predicted by the dotted expectation line 
in Fig. 2 after the introduction of an adiabatic term [19] 
to account for the kinematic cutoff. 

In trying to find a solution to these wee conundrums, 
this author searched [16, 17] for a description of the 
method of alignment between 2.8 MeV neutron source, 
LXe cell, and scintillator cell used to detect scattered 
neutrons (the scattering angle defining the recoil energy 
in the measurements). None was found, and also an ab- 
sence of treatment for the uncertainty in recoil energy 
that any angular misalignment would bring, therefore un- 
realistically assumed to be nil 8 . Based on the geometry 
of this experiment, size of the detector cells and the rel- 
atively short distance between them, a modest misalign- 
ment in the relative position of the cells and source could 
enlarge the recoil energy error bars significantly and/or 
shift the datapoints in Fig. 2 in energy: as can be seen 
from the inset in Fig. 2, a cumulative misalignment by as 
little as 5°, which amounts in the experiment by Manzur 
et al. to 2-3 cm of cumulative displacement in the as- 
sumed relative position between the detector cells and/or 
between them and the neutron source, is sufficient to re- 
move the paradoxical situations described above (e.g., 
the lowest energy datapoint in Fig. 2 shifts from 4 keV r 
to 6 keV r over 5°). A certain discomfort arises from 
the simultaneous realization that the XENON10 light- 
WIMP sensitivity claimed in [28] using the £ e ff measure- 
ments from Manzur et al. is also critically dependent 
on this matter: the same unaccounted-for displacements 
by few cm would result, quite literally, in a decrease in 
XENON10 light- WIMP sensitivity by orders of magni- 
tude. While this should be reason enough to take the 
limits claimed in [28] with caution, much graver concerns 
about the methodology used by Manzur et al. and earlier 
authors to arrive to £ e ff and Q y will be expressed in Sec. 
III. Considerably less information (trigger and software 
efficiency, S2 gain, analysis procedure, etc.) is provided 



As it turns out, no goniometric equipment was employed by 
Manzur et al. The approximate scattering angles were derived 
from photographs of the setup taken from the vertical [27]. In 
the opinion of this author, this method does not guarantee the 
desired angular accuracy discussed in the main text. 
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by Manzur et al. in [16, 17] on the subject of their deriva- 
tion of Q y , as compared to their very extensive treatment 
of C e s , preventing an additional investigation here of how 
the neglected systematic effect described in Sec. Ill may 
have affected the Q y datapoints in Fig. 2. A skewness 
similar to that affecting the low-energy SI light peaks 
discussed in Sec. Ill can be noticed in the single S2 light 
example shown by Manzur et al. (Fig. 8c in [16]), possi- 
bly indicating the presence of S2 threshold effects at the 
level of Q ~few electrons. The reader is referred to Sec. 
Ill for more details. 

As a last remark on the difficulties in the determination 
of light- WIMP sensitivity via S2 light, the bottom panel 
in Fig. 1 shows the effect of adopting the expected en- 
ergy scale (central value of dotted lines in Fig. 2) rather 
than that generated by the latest version of the best- 
fit Monte Carlo, taking into account the changes in bin 
width that the energy translation results in. It would be 
possible to obtain a near-perfect agreement between the 
example light WIMPs in the figure (intentionally cho- 
sen in the region of interest for DAMA, CoGeNT and 
CRESST [6]) and this spectrum, with the introduction 
of the adiabatic term to account for the expected kine- 
matic cutoff or by examining other uncertainties in the 
expectations. It is sobering to realize that one can go 
from this bizarre coincidence to severe exclusion of the 
candidates over a shift in a notoriously ill-defined en- 
ergy scale by a mere ~ 3 keV r . While the presence of a 
high-rate, homogeneously-distributed signal at very low- 
energy in the XENON 10 or XENON 100 detector would 
be hard to explain in an efficient self-shielding medium 
like LXe without invoking weakly-interacting particles, 
any such possible agreement with these other detectors 
should be presently de-emphasized: at the time of this 
writing, any low energy scale for LXe should be consid- 
ered highly speculative, such is the paucity of reliable 
experimental data and detailed models of recoil response 
available for LXe. 



Ill: UNCERTAINTY IN QUENCHING FACTOR 

The information presently available on scintillation 
and ionization yield from low-energy nuclear recoils in 
LXe is plagued by contradictory observations and inter- 
pretations. In particular, £ c g is expected to decrease 
with decreasing recoil energy [8] for the reasons given 
above, but published measurements display all possi- 
ble behaviors. Recent attempts to reconcile these [29] 
fall short, lacking in a critical examination of the anal- 
ysis methods and experimental techniques utilized in 
the existing calibrations. In order to elucidate some of 
these issues, extensive MCNP-PoliMi [30] simulations of 
the most recent assessment of C e fi> by Manzur et al. 
in [16, 17], have been performed here. This examina- 
tion has revealed several issues affecting this type of 




FIG. 3: MCNP-PoliMi simulated distribution of recoil en- 
ergies from different components of the spectrum measured 
by Manzur et al, comparable to Fig. 9 (top) in [16]. The log 
scale generates subtle differences in appearance between these 
figures. 



measurement 9 . Attention was paid to follow the pre- 
scriptions given in [16, 17] regarding energy resolution, 
SI light yield at different drift field values, cuts affect- 
ing the data, effect of threshold efficiency, etc., obtain- 
ing a good agreement with the particulars provided in 
[16, 17] on distribution of recoil energies, contributions 
from different types of events (inelastic scatters, multiple 
scatters, scattering on inert materials, etc.), and neutron 
time-of-flight (TOF) information. In this respect, Fig. 3 
here allows a direct comparison to Fig. 9 in [16]. 

An important limitation to the measurements by 
Manzur et al. is not readily appreciated from a cursory 
inspection of [16]: the position of the peak or maximum 
in the measured distributions of SI scintillation light for 
neutron scattering angles below 6 ~ 60° is entirely de- 
fined by the effect of the threshold efficiency (software 
and triggering). In other words, for recoil energies below 
^20 keV r , the left shoulders responsible for the notice- 
able skewness in distributions of SI light like those in 
Fig. 8a in [16] and Fig. 6.3 in [17] are solely the result of 
this threshold efficiency: the true scintillation maximum 
is expected at values much lower than the position of 
this peak. These true maxima are not reachable due to 
an scintillation yield (4.3 to 10.8 photoelectrons per keV 
of electron equivalent energy, PE/keV ee , depending on 
drift field and mode of operation) that is insufficient for 



Several of these issues are minor, and yet illustrative of the lack 
of scrutiny with which some of these results have been embraced. 
For instance, simple geometrical considerations suffice to demon- 
strate that the distance between LXe and scintillator cells quoted 
in [16] (16-20 cm, information not mentioned in [17]) would lead 
to an uncertainty in the recoil energies probed almost three times 
larger than claimed. This was confirmed by initial simulations. 
This issue is traceable to a mistake in units of distance in [16], 
where inches should have been used [27]. 
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FIG. 4: Top: Simulated distributions of SI light in the exper- 
imental setup of Manzur et al. [16, 17], for 5.7 keV r recoils. 
Black histograms use the nominal C e g claimed in [16, 17] for 
this recoil energy, red histograms a logarithmic fit to their 
datapoints above 20 keVr (dashed line in Fig. 5 bottom, see 
text). The effect of the threshold efficiency is to shift the 
peak or maximum in the distributions to an artificially larger 
value. Both £ c fj hypothesis remain distinguishable through 
their predicted event rate (area under solid histograms). Bot- 
tom: The ability to distinguish these two largely different 
values of £ e g is lost once a normalization to the experimental 
rate (data points from Fig. 8a in [16]) is performed, an un- 
fortunate method of analysis followed by Manzur et al. and 
earlier authors [15, 31, 32]. 




10 100 
recoil energy (ke\^ ) 



FIG. 5: Top: Position of the peak or maximum in SI scintilla- 
tion (units of electron equivalent energy) observed by Manzur 
et al. as a function of recoil energy [33]. The dashed line is 
derived from the logarithmic fit in the bottom panel, under 
the premise that the ratio of peak keV ee to keVr should be 
equal to C c s in ideal experimental conditions (see text). Bot- 
tom: C e g datapoints from Manzur et al. together with the 
yCeff derived by the ZEPLIN collaboration (red band [34]). 
The dashed line is a logarithmic fit to datapoints not affected 
by threshold efficiency (> 20keV r ). Lower energy datapoints 
deviate from the fit by the effect of the threshold efficiency 
(see text). Cutoffs to this fit anywhere in the range 3-6 keV r 
produce fits to the experimental data similar in quality to 
those shown in the bottom panel of Fig. 4. 



an optimal exploration of such low recoil energies. For 
instance, assuming a value of £ c ff =0.1 at 6 keV,-, this 
peak would then be expected at between 2.6 and 6.5 PE, 
whereas the effect of the threshold efficiency starts to be 
noticeable already at ~15 PE for both the single and 
two-phase modes of operation. The threshold efficiency 
is provided in Fig. 10 in [16] and top of Fig. 4 here. 

This is not an insurmountable limitation, as long 
as careful attention is paid to the comparison between 
Monte Carlo simulations and data that ultimately leads 
to a best-fit £ e ff for each scattering angle. Unfortunately, 
in the experiments described by Manzur et al. no at- 
tempt was made to include a comparison between the 
expected recoil rate and that observed 10 . Instead, the 
simulated distributions of SI light yield were normalized 



to the observed event rate prior to the chi-square analysis 
leading to the best-fit value for £ gf, a step mentioned in 
[16] and described with more detail in Sec. 6. 1. 2 of [17]. 

The nature of the problem created by this form of data 
analysis is illustrated in Fig. 4, using the 6 keV r SI single- 
phase measurement provided as an example in [16]. The 
top panel displays the expected SI light yield distribution 
for the nominal £ e ff =0.078 obtained by Manzur et al. at 
this energy 11 , derived from the present simulation follow- 
ing the data cuts and energy resolution prescribed in [17], 
before (dotted histograms) and after (solid histograms) 
including the effect of threshold efficiency. Its equivalent, 
but using a considerably smaller and energy-dependent 
£ e g derived from a logarithmic fit (Fig. 5, bottom) to 
£ e ff datapoints with E rec >20 keV r (i.e., those not af- 



This was prevented by difficulties in controlling the stability of 
the DD neutron generator output towards the end of its target 
life, and by throughput limitations of the data acquisition system 
during runs at the smallest scattering angles [27]. 



£ e fT is expected to be an energy-dependent quantity. Use of a 
constant value across the measured SI spectrum during the chi- 
square analysis, as done in [16, 17] is less than ideal. An effort 
was nonetheless made in [17] to examine the effect of a linear 
energy dependence around the best-fit C e g. 
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fected by the effect of threshold efficiency) is also dis- 
played. As can be observed from the solid histograms in 
the top panel of Fig. 4, these two possibilities lead to dis- 
tinguishable predictions, specifically a decrease in recoil 
rate per unit exposure to the neutron source by 40% in 
the second scenario. However, once the normalization to 
the experimental rate is performed (Fig. 4, bottom), this 
information is lost and both £ e g scenarios yield a com- 
parable good fit to the experimental data. In conclusion, 
the method of analysis followed in [16, 17] (and earlier 
experiments discussed next) cannot distinguish between 
these two largely different values of £ e g. 

Useful additional information was presented in an un- 
published progress report [33] and reproduced here in 
the top panel of Fig. 5, where the position of the ob- 
served peak or maximum in SI light yield (in units of 
keV ee rather than PE) is displayed as a function of mean 
recoil energy for early measurements by Manzur et al. 
For data points not affected by the threshold efficiency 
(>20 keV r ) the ratio between both energies is, as ex- 
pected by the definition of C c s, close to the value of 
£ e s eventually obtained (Fig. 5, bottom). For smaller 
recoil energies, for which the threshold efficiency affects 
and defines the position of the peak, this approximation 
breaks down, as a result of the true position of this peak 
having been artificially "pushed up" in energy, the ef- 
fect illustrated here in Fig. 4 (top). This can be ob- 
served in, for instance, Figs. 11 a,c in [16], where the 
chi-square minimum appears at £ e g ^0.055, rather than 
at the larger ~0.12 that would be naively expected from 
said ratio and a 10.8 PE/keV ee light yield. This trend 
in [16, 17] for the chi-square minimum to appear well- 
below an expectation based on the (artificial) position 
of the scintillation peak, is in itself pointing at a rapidly 
decreasing £ e g towards zero recoil energy. But this is rec- 
ognized in [16, 17] and not the point at stake: instead, 
the top panel in Fig. 5 reveals that the measurements 
by Manzur et al. not affected by the threshold efficiency 
may very well be pointing at a vanishing £ e ff somewhere 
in the few keV r region, a behavior not dissimilar to that 
claimed by the ZEPLIN collaboration [34]. The change 
in slope in Fig. 5 (top) appears exactly where expected 
from the limitations imposed by threshold efficiency and 
light yield in this experiment 12 . It is inevitable to sus- 
pect that the combination of insufficient light yield and 
regrettable normalization of simulated to experimental 
rates used by Manzur et al. and others [15, 31, 32] may 
have biased £ c g towards considerably larger values than 



12 The asymmetry in the SI distribution due to the effect of thresh- 
old efficiency is already clearly visible in the 10 keV r data in Fig. 
6.3 of [17]. At 20 keV r and assuming a C c g ~0.1-0.15, at least 
the lower light yield, non-zero drift field measurements are ex- 
pected to be affected (the C e g datapoints of Manzur et al. are 
an average of all conditions tested for each angle). 
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FIG. 6: Deficit in measured vs. expected single Cs,I recoil 
rate under irradiation of a CsI[Na] scintillator with 2.8 MeV 
monochromatic neutrons. Vertical arrows mark the value of 
the kinematic cutoff in recoil energy expected for ionization 
and scintillation in this medium ([8], see text). 



those to be measured in more ideal conditions, or at least 
when properly taking into account systematic differences 
between expected and observed signal rates. The loga- 
rithmic fit to the points above 20 keV r employed to illus- 
trate this discussion (Fig. 5, bottom) lays slightly beyond 
the one-sigma uncertainty in £ e g quoted by Manzur et 
al. 13 The adoption of a similar £ e g would dramatically re- 
lax claimed constraints on light- WIMPs from XENON10 
and XENON100 [28, 35]. 

An attempt has not been made here to simulate earlier 
measurements of £ e g as a function of neutron scattering 
angle. However, at least in the analysis of the experi- 
ments described in [15, 31, 32], not only the overall nor- 
malization in rate was left as a free parameter (enabling 
the same deleterious effects illustrated here by Fig. 4), 
but additional degrees of freedom were introduced into 
the fits by way of a background model intended to ac- 
count for neutron scattering in inert components, largely 
dominant in those other geometries 14 . This is clearly 



13 An unquantified very large increase in signal rate towards small 
scattering angles was noticed in the measurements by Manzur et 
al. [27]. Present simulations predict that this rate should have 
increased by a factor ~200 in going from =125° (67 keV r ) to 
9 =25° (4 keV r ), before accounting for threshold effects. Yet 
only a small deficit (~40%, Fig. 4 top panel) in rate is expected 
from an £ e ff much smaller than that determined by Manzur et 
al, leaving ample margin for a severe misidentification of the 
low-energy C e g, even more dramatic than what has been con- 
templated here. 

14 The detector used by Manzur et al. is properly designed for the 
job at hand, with only a bare minimum of inert parts around the 
LXe cell. The largely dominant "materials background" in [31, 
32] arises from a suboptimal detector design, causing neutrons to 
have a very small probability of reaching and exiting the active 
cell without interacting on inert materials, creating a background 
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a flawed approach to the analysis of such experimental 
data. The dismayingly poor quality of the calibration 
data in [31, 32] is worth emphasizing at this point: an in- 
terested reader is invited to inspect the discussion around 
Fig. 4-10 in [32] for crucial steps in the data analysis ob- 
viated in [31]. In the opinion of this author and others 
[27], once the obvious residual threshold effect surviving 
the subtraction of accidentals is corrected for, there is 
little to no evidence in the dataset of [31, 32] for any us- 
able calibration information at and below 10 keV r . This 
calibration played a central role in the attempted justifi- 
cation of recent XENON100 limits [35]. 

To finalize, Fig. 6 is offered to illustrate the reality of 
the concerns expressed in this paper. It displays the rate 
of single Cs or I recoils in a CsI[Na] scintillator irradiated 
by a 2.8 MeV DD neutron generator in a setup similar to 
that used by Manzur et al. These data are part of an on- 
going series of precision quenching factor calibrations for 
CsI[Na] , CsI[Tl] and NaI[Tl] at the University of Chicago, 
employing combined electron recoil (Compton scattering) 
and nuclear recoil measurements on a single goniometric 
table. While these recent data should be considered pre- 
liminary, a large deficit in signal rate is observed, with 
an onset coincident with the predicted kinematic cutoff 
for Cs,I recoils in CsI[Na] [8]. At the time of this writ- 
ing, data acquisition triggering throughput and threshold 
efficiency effects are found to be far from sufficient to ac- 
count for its magnitude (the measured low-energy light 
yield of this scintillator in this setup is —9 PE/keV ee ). 
Following the incorrect analysis procedure criticized in 
this section, a monotonically increasing quenching factor 
towards zero recoil energy would be naively derived from 
the CsI[Na] data at hand. The lowest energy data point 
centered at 5 keV r would be assigned a quenching factor 
of —20%, with the highest energy points in agreement 
with earlier measurements [36] at ~8%. With a correc- 
tion in place for this systematic deficit in rate, the derived 
quenching factor should be much smaller for the lowest 
energies measured. An upcoming publication will discuss 
the implications of these measurements for the exact lo- 
cation of the DAMA favored region in light- WIMP phase 
space [6, 8], impact on efforts to reproduce the annual 
modulation effect using CsI[Tl] crystal arrays [37], and 
prospects for a low-energy neutrino measurement using 
CsI[Na] at an intense neutron spallation source [38]. 



IV: CONCLUSIONS 

Even after considering the difficulties described in the 
introduction, common to all techniques, it is possible to 
distinguish between detector technologies more adapted 



that completely swamps the signals sought. 



to a search for light WIMPs than others, based exclu- 
sively on the existing knowledge of the response of the 
detector in the few keV recoil energy region. Two fami- 
lies of devices can be examined: those for which a reliable 
low-energy scale exists, and those for which it is presently 
lacking. Belonging to the first group, DAMA/LIBRA 
Nal(Tl) scintillators benefit from a convenient low-energy 
signal at 3.2 keV originating in a 40 K internal contam- 
ination, able to anchor the energy scale near threshold, 
even if a relatively small uncertainty still remains in the 
quenching factors that define the recoil energy scale [6, 8]. 
CDMS germanium bolometers [39] profit from peaks at 
1.3 keV and 10.4 keV from 71 Ge activation following pe- 
riodic neutron calibrations. The quenching factor mea- 
sured in CDMS-Ge is similar to theoretical expectations 
and the observations from independent germanium ex- 
periments. CRESST bolometers are able to use narrow 
lines of known origin at 3.6 keV and 8 keV for a low- 
energy reference, having measured the quenching factors 
for the recoiling species in their crystals using a num- 
ber of techniques [5]. CoGeNT detectors benefit from a 
number of known narrow cosmogenic lines in the range 
1.1-11.1 keV and an optimal linearity and energy resolu- 
tion. Their quenching factor has been measured in a ded- 
icated reactor experiment down to sub-keV recoil ener- 
gies, obtaining an excellent agreement with expectations 
and other measurements [40]. Belonging to the second 
group we can identify CDMS silicon bolometers: spatial 
corrections are known to affect their recoil energy scale, 
generating a large disagreement with theoretical values 
of the quenching factor (even larger when compared to 
other experimental values), and a potentially large shift 
in recoil energy [6]. At the time of this writing, detectors 
based on LXe lag behind all other dark matter detection 
techniques in a knowledge of the energy scale and of the 
mechanisms governing signal production by few keV r nu- 
clear recoils. Recent efforts to provide an electron-recoil 
low-energy calibration reference via a 83m Kr source [41] 
are a welcome step in the right direction. 

To summarize, important flaws in the methodology 
used to determine the value of £ e ff, Qy and the recoil en- 
ergy scale in liquid xenon detectors have been described. 
Data-quality concerns affecting recent calibrations, such 
as those in [31, 32], have been pointed out. Once these 
issues are addressed, an agreement with the expected de- 
crease in ionization and scintillation yield from nuclear 
recoils below kinematic cutoff should be obtained for this 
medium. In particular, the traditionally observed deficit 
in low-energy response to AmBe neutron-induced recoils 
in LXe should be easier to understand by the inclusion 
of well-known radiation effects presently being ignored. 
A proper control of the systematics affecting monochro- 
matic neutron scattering measurements of C e g and Q y 
has been shown to be lacking: bold as this may sound, 
none of the measurements of this type performed up until 
now below -10 keV r by the XENON10 or XENON100 
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collaborations can be assigned much worth. Separately, 
an effort should be made to evolve the "best-fit Monte 
Carlo" method to include radiation effects known to me- 
diate the generation of signal carriers (free electrons, di- 
rect scintillation) before an attempt to match simulations 
and data is made. Proper account of the uncertainty in 
the recoil energy scale should be taken before attempting 
to extract limits from spectra generated by this method. 
In its present naive form, it is not possible to defend it. 

Markedly different values of C e g and Q y should be 
expected from improved LXe calibration methodologies. 
The recoil energy scale of ionization (S2 light) spectra 
should also be affected by the proposed improvements, 
as described in the commentary around Figs. 1 and 2. 

In view of the numerous issues raised here, one is forced 
to conclude that recent attempts to extract light- WIMP 
sensitivity from XENON10 and XENON100 data [28, 35] 
arc premature and overly optimistic. The sensitivity 
of present LXe detectors to WIMPs heavier than ~10 
GeV/c 2 is also affected by the considerations presented 
here, albeit to a lesser (but calculable) extent. 

The author is indebted to D.N. McKinsey for many 
frank exchanges that made the present analysis possible, 
and for calling his attention to the limitations affecting 
the measurements in [31, 32]. D. Hooper, N. Weiner and 
K. Zurek provided useful suggestions. 
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